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Complex dynamical networks consisting of many components that 
interact and produce each other are difficult to understand, espe- 
cially, when new components may appear. In this paper we outline 
a theory to deal with such systems. The theory consists of two parts. 
The first part introduces the concept of a chemical organization as 
a closed and mass-maintaining set of components. This concept al- 
lows to map a complex (reaction) network to the set of organizations, 
providing a new view on the system's structure. The second part 
connects dynamics with the set of organizations, which allows to 
map a movement of the system in state space to a movement in the 
set of organizations. 

Our world is changing, qualitatively and quantitatively. The characteristics 
of its dynamics can be as simple as in the case of a friction-less swinging 
pendulum, or as complex as the dynamical process that results in the creative 
apparition of novel ideas or entities. We might characterize the nature of a 
dynamical process according to its level of novelty production. For example, 
the friction-less swinging pendulum implies a process where the novelty is only 
quantitative. Whereas the process of biological evolution is highly creative and 
generates qualitative novelties, which then spread in a quantitative way. 

Fontana and Buss [T| called processes and systems that display the produc- 
tion of novelty, constructive (dynamical) processes and constructive (dynami- 
cal) systems, respectively. Constructive systems can be found on all levels of 
scientific abstraction: in nuclear physics, where the collision of atoms or sub- 
atomic particles leads to the creation of new particles; in molecular chemistry, 
where molecules can react to form new molecules; or in social systems, where 
communication can lead to new communication [2]. As a result of a combina- 
torial explosion, it is easy to create something that is new, e.g., a molecule or 
a poem that is unique in the whole known universe. Nevertheless, it should be 
noted that novelty is relative. Whether something is considered to be new or 
not, depends on what is already there. Thence it follows that whether a system 
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appears as a constructive dynamical system is usually context dependent. 

Despite the fact that a large amount of interesting dynamical processes are 
constructive, classical systems theory does not, conveniently, take novelty into 
account. Classical systems theory assumes a given static set of components. 
For example, a classical systems analysis would first identify all components 
and their relations; then would identify a state space with fixed dimensionality, 
e.g., the state of our pendulum could be described by a two-dimensional vector 
specifying angle and angular velocity. 

The lack of a theory for constructive dynamical systems has been identi- 
fied and discussed in detail by Fontana and Buss P in the context of a theory 
for biological organization. As a solution, they suggested the important con- 
cept of a (biological) organization as an operationally closed and dynamically 
self-maintaining system. Taking this idea further, we suggested to divide the 
theoretical approach in a static and a dynamic analysis [3]. Here, we present 
the results of our study for general reactions systems, integrating stoichiom- 
etry and the concept of mass-maintenance in the definition of organizations, 
which is a prerequisite for a broad applicability. 

Reaction Systems. The theory described herein aims at understanding re- 
action systems. A reaction system consists of reactants, which we refer to as 
molecules. An interaction among molecules that lead to the appearance or 
disappearance of molecules is called a reaction. In particular, reaction systems 
are used to model chemical processes, but in general, they can be applied to 
virtually any domain where elements interact to influence the production of 
other elements, e.g., population dynamics, evolution of language, economy, or 
social dynamics. 

We have to distinguish between a reaction system as an abstract descrip- 
tion of all possible molecules (and their reactions), and a reaction vessel, which 
contains concrete instances of molecules from the set of all possible molecules. 
In general, the description of a reaction system can be subdivided into three 
parts: (1) the set of all possible molecules A4, (2) the set of all possible 
reactions among all the possible molecules TZ, and (3) the dynamics, which 
describes how the reactions are applied to a collection of molecules inside a 
reaction vessel. 

Static Analysis 

Algebraic Chemistry. In the first part of the paper, we are only concerned 
with the static structure of a reaction system, that is, the molecules and the 
reactions. And instead of considering a state (e.g., a concentration vector), we 
limit ourself to the analysis of the set of molecules present in that state. We 
introduce the concept of an algebraic chemistry, which is a reaction network, 
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including stoichiometric information, from which we will derive the organiza- 
tional structure of the system. 

Definition 1 Given a set A4 of elements (called molecules) and a set of reaction 
rules given by the relation TZ : Vm{M) x Vm{M). We call the pair (M,TZ) 
an algebraic chemistry. 

Vm(C) is the set of all multisets with elements from C. A multiset differs 
from a set in the fact that the same element can appear more than one time. 
The frequency of occurrence of an element a in a multiset A is denoted by 
#(a G A). For simplicity, we adopt a notion from chemistry to write reaction 
rules. Instead of writing ({s±, s 2 , ■ ■ ■ , s n }, {s[, s' 2 , . . . , s' n ,}) G TZ we write: Si + 
s 2 + ■ ■ ■ + s n — > s[ + s' 2 + • ■ ■ + s' n , . Given the left hand side molecules 
A = {si, s 2 , ■ ■ ■ , s n } and the right hand side molecules B = {s[, s' 2 , ■ ■ ■ , s' n ,}, 
we write (A — > B) G TZ instead of (A, B) G TZ. A — > B represents a chemical 
reaction equation where A is the multiset of molecules on the left hand side 
(also called reactants) and B the multiset of molecules on the right hand side 
(also called products). 

Input and Output. There are many processes that give rise to an inflow and 
outflow, such as, incident sunlight, decaying molecules, or a general dilution 
flow. In this paper we handle input simply by adding a reaction rule — > a 
to the reaction rules TZ for every molecule a that is in the inflow. denotes 
the empty set. Equivalently, for an output molecule a (e.g., a molecule that is 
decaying) we add the rule a — > to the reaction rules TZ. 

Semi-Organization. In classical analysis, we study the movement of the 
system in state space. Instead, here, we consider the movement from one 
set of molecules to another. As in the classical analysis of the dynamic of 
the system, where fixed points and attractors are considered more important 
than other states, some sets of molecules are more important than others. 
In order to find those sets, we introduce some properties that define them, 
namely: closure, self-maintenance, semi-organization, mass-maintenance, and 
finally the organization. 

All definitions herein refer to an algebraic chemistry (At, TZ). The first re- 
quirement, called closure, assures that a set of molecules contains all molecules 
that can be produced by reactions among those molecules. 
Definition 2 A set C C M. is closed, if for all reactions (A — > B) G TZ, with A 
a multiset of elements in C (A G Vm{C)), B is also a multiset of elements in 
C (Be V M (C)). 

Given a set S C M., we can always generate its closure Gcl(S) according 
to the following definition: 

Definition 3 Given a set of molecules S C M, we define Gcl(S) as the smallest 
closed set C containing S. We say that S generates the closed set C = Gcl(S) 
and we call C the closure of S. 



3 



The closure implies a union and an intersection operator on the closed sets 
Oql of an algebraic chemistry. Given two closed sets U and V , the closed set 
generated by their union (U Ucl V) and intersection (U \1cl V) is defined as: 
UU CL V = G C l{U U V), and U VA CL V = G CL (U n V), respectively. Triv- 
ially, closed sets form a lattice (0 C l, Ucl, ^cl), which is a common algebraic 
structure (a poset in which any two elements have a greatest lower bound and 
a least upper bound, see supl. mat.). The property of closure is important, 
because the closed set represents the largest possible set that can be reached 
from a given set of molecules. Furthermore a set that is closed cannot generate 
new molecules and is in that sense more stable than a set that is not closed. 

The next important property, called self-maintenance, assures that every 
molecule that is used-up within a set, is produced within that set. 
Definition 4 A set of molecules S C Ai is called self-maintaining, if all molecules 
s G S that are used-up within S are also produced within that set S. 
More precisely, we say that a molecule k G M. is produced within a set C C M., 
if there exists a reaction (A -> B), with A G V M {C), and G A) < #(& G 
5). In the same way, we say that a molecule k G C is used-up within the set 
C, if, within the set C, there is a reaction (A — > £?) with A G Vm(C), and 

#^)>#e5). 

Taking closure and self-maintenance together, we arrive at the important 
concept of a semi-organization. 

Definition 5 A semi-organization O C .M is a set of molecules that is closed 
and self-maintaining. 

Organization. In a semi-organization, all molecules that are used-up are 
produced; but, this does not guarantee that the total amount of mass can be 
maintained. A small, but important example is the reversible reaction in a 
flow reactor: M = {a,b},1Z = {a — > b,b — > a, a — > 0,6 — > 0}, where both 
molecules also decay. O = {a, b} is a semi-organization, because the set is 
closed, a is produced by the reaction b — > a, and 6 is produced by the reaction 
a — > 6. But, obviously, the system {a, 6} is not stable, because both molecules 
decay and are not sufficiently reproduced, so that they will finally vanish. The 
solution to this problem is to consider the overall ability of a set to maintain 
its total mass. We call such sets mass-maintaining: 

Definition 6 Given an algebraic chemistry (At,lZ) with m = \A4\ molecules 
and n = \1Z\ reactions, and let M = (m^j) be the (m x n) stoichiometric 
matrix implied by the reaction rules TZ, where rriij denotes the number of 
molecules of type % produced in reaction j. A set of molecules C C M is called 
mass-maintaining, if there exists a flux vector v e 1" such that the three 
following condition apply: (1) for all reactions (A — > B) with A G Vm{C) the 
flux V(a->b) > 0; (2) for all reactions (A — > B) with A Vm{C), v^a^b) = 0; 
and (3) for all molecules % G C, fi > with (fi, . . . , f m ) = Mv. 
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V(a-^b) denotes the element of v describing the flux of reaction A — > B. For 
the example above, the stoichiometric matrix becomes M = ((—1, 1), (1, — 1), 
(—1,0), (0, —1)), and we can see that there is no positive flux vector v e M 4 , 
such that Mv > 0. In fact, only the empty semi-organization {} is mass- 
maintaining. In case a and b would not decay, 1Z = {a — > b, b — > a}, the set 
{a, b} would be (as desired) mass-mainatining, because threre is a flux vector, 
e.g., v = (1.0, 1.0), such that Mv = > with M = ((-1, 1), (1, -1)). Now, 
closure and mass-maintenance together lead to the central definition of this 
work: 

Definition 7 A set of molecules O C M. that is closed and mass- maintaining 
is called an organization. 

An organization represents an important combination of molecular species, 
which are likely to be observed in a reaction vessel on the long run. A set 
of molecules that is not closed or not mass-maintaining would not exists for 
long, because new molecules can appear or some molecules would vanish, re- 
spectively. The condition "mass-maintaining" is stronger and more difficult to 
compute than the condition "self-maintaining". In fact, the former implies the 
latter, because in a mass-maintaing set, a molecule % that is used-up within that 
set must also be produced within that set in order to achieve a non-negative 
production rate /j. Thus we can say: 
Lemma 1 Every organization is a semi-organization. 

Example 1 In order to illustrate the new concepts, we consider a small exam- 
ple where there are just four molecular species Ai = {a,b,c,d}, which react 
according to the following reaction rules 1Z = {a + b — > a + 2b, a + d — > 
a + 2d, b + c — > 2c, c — > b, b + d — > c, b — > 0, c — > 0, d — > 0}, where b, c, d decay 
spontaneously and a is a permanent catalyst. Although the reaction system is 
small, its organizational structure is already difficult to see when looking at the 
rules or their graphical representation (Fig. 1). In Fig. 1, all 16 possible sets 
of molecules are shown as a lattice. There are 9 closed sets, 8 self-maintaining 
sets, 7 mass-maintaining sets, and 7 semi-organizations, 6 of which are orga- 
nizations. The organizations are the only combination of molecules that can 
reside in a reaction vessel for a long time. We can also see immediately that 
a reactor containing {b, c, d} must have a transient dynamics "down" where 
molecular species are lost, and that a reactor containing {a, b, d} will have a 
transient where a new molecular species will appear, and so on... 

Different reaction systems 

Finding all organizations of a general reaction system appears to be compu- 
tationally difficult. One approach is to find the semi-organizations first, and 
then check, which of them are also mass-maintaining. The property of the 
set of organizations and semi-organizations depends strongly on the type of 
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system studied. We will present here four types of systems and analyze their 
properties with respect to the newly introduced concepts. 

Catalytic Flow System . In a catalytic flow system all molecules are used- 
up by first-order reactions of the form {k} — » (dilution) and there is no 
molecule used-up by any other reaction. So, each molecule k decays sponta- 
neously, or equivalently, is removed by a dilution flow. Apart from this, each 
molecule can appear only as a catalyst (without being used-up). Examples 
of catalytic flow system are the replicator equation jl], the hypercycle the 
more general catalytic network equation [6 , or AlChemy [7j. Furthermore some 
models of genetic regulatory networks and social system jHj are catalytic flow 
systems . 

Lemma 2 In a catalytic flow system , every semin-organization is an organiza- 
tion. 

In a catalytic flow system we can easily check, whether a set O is an organi- 
zation by just checking whether it is closed and whether each molecule in that 
set is produced by that set. Furthermore, given a set A, we can always gen- 
erate an organization by adding all molecules produced by A until A is closed 
and than removing molecules that are not produced until A is self-maintaining. 
With respect to the intersection and union of (semi-)organizations the set of 
all (semi-) organizations of a catalytic flow system forms an algebraic lattice 
(see below), which has already been noted by Fontana and Buss PQ. 

Reactive Flow System . In a reactive flow system all molecules are used- 
up by first-order reactions of the form {k} — > (dilution). But as opposed 
to the previous system, we allow arbitrary additional reactions in 1Z. This 
is a typical situation for chemical flow reactors or bacteria that grow jU]. In 
a reactive flow system , semi-organizations are not necessarily organizations, 
as in the reversible reaction example shown before. Nevertheless, both, the 
semi-organizations and the organizations form a lattice (O, U,n), Moreover, 
the union (U) and intersection (n) of any two organizations is an organization 
(see below). 

Consistent Reaction System . In a consistent reaction system there are 
two types of molecules: persistent molecules P and non-persistent molecules. 
All non-persistent molecules k e M \ P are used-up (as in the two systems 
before) by first-order reactions of the form {k} — > B with k ^ B\ whereas a 
persistent molecule p G P is not used-up by any reaction at all. An example 
of a consistent reaction system is Example 2, where a is a persistent molecule. 
The consistent reaction system is the most general of the four systems where 
the semi-organizations and organizations always form a lattice, and where the 
generate organization operator can properly be defined (see below). As in a 
reactive flow system , not all semi-organizations are organizations. 
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General Reaction System . A general reaction system consists of arbi- 
trary reactions. Examples are planetary atmosphere chemistries JU]- m a 
general reaction system , neither the set of organization nor the set of semi- 
organizations necessarily form a lattice. Because of the lack of properties there 
is an actual difficulty in studying those systems. The analysis gets much easier, 
if we are able to transform the system into something else, e.g., a consistent 
reaction system . This could be done by introducing a small outflow for each 
molecule. 

Common Properties of Consistent Reaction System Consistent reac- 
tion systems (including reactive flow system and catalytic flow system ) possess 
some comfortable properties that allow us to present a series of useful defini- 
tions and lemmas. In a consistent reaction system , given a set of molecules C, 
we can uniquely generate a self-maintaining set, a semi-organization, a mass- 
maintaining set, and an organization in a similar way as we have generated 
a closed set. And like for closed sets, we can define the union and intersec- 
tion on self-maintaining sets, semi-organizations, mass-maintaining sets, and 
organizations, respectively. Furthermore, each, the self-maintaining sets, semi- 
organizations, mass-maintaining sets, and organizations form a lattice together 
with their respective union and intersection operators. This does not gener- 
alize to general reaction system , because for a general reaction system we 
cannot uniquely generate a self-maintaining set, a semi-organization, a mass- 
maintaining set, nor an organization as in the case of a consistent reaction 
system . 

Definition 8 Given a set of molecules C C M, we define G$m(C) as the biggest 
self-maintaining set S contained in C. We say that C generates the self- 
maintaining set S = Gsm{C). 

In order to calculate the self-maintaining set generated by C, we remove those 
molecules that are used-up and not produced within C, until all molecules used- 
up are also produced, and thus reaching a self-maintaining set. The operator 
Gsm (generate self-maintaining set) implies the union Usm and intersection 
I~l5jv/ on self- maintaining sets: Given two self maintaining sets S± and S 2 , the 
self- maintaining sets generated by their union (Si U C l S 2 ) and intersection 
(Si T\ C L S 2 ) are defined as: Si U 5A/ S 2 = G SM (Si U S 2 ), and Si n SM S 2 = 
GsM(Sir\S2), respectively. And as already mentioned, in a consistent reaction 
system , (Osm,^sm,^sm) forms a lattice, where Osm is the set of all self- 
maintaining sets of an algebraic chemistry. Finally note that, if S is self- 
maintaining, its closure Gcl(S) is self-maintaining, too (in consistent reaction 
systems ). 

There are many ways in which we can generate a semi-organization from 
a set. We will present here the simplest one, which implicitly assumes that 
molecules are produced quickly and vanish slowly. This assumption leads to 
the largest possible semi-organization generated by a set: 
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Definition 9 Given a set of molecules C C M., we define GsoiC) as Gsm(Gcl(S)). 
We say that C generates the semi- organization O = GsoiC). 

In the same way as before, the generate semi-organization operator Gso{C) 
implies the union Uso and intersection Uso on semi-organizations, namely 
1 U so 2 = G SO (Oi U 2 ), and 1 n so 2 = G SO (Oi n 2 ), respectively, 
which implies the lattice of semi-organizations. 

Definition 10 Given a set of molecules C C M. ) we define Gmm(C) as the 
biggest mass-maintaining set S contained in C. We say that C generates the 
mass-maintaining set S = Gmm{C). 

For consistent reaction systems , Gmm(C) is always defined, because the union 
(U) of two mass-maintaining sets is mass-maintaining; and further, every set is 
either mass-maintaining, or it contains a unique biggest mass-maintaining set. 
Thus from every set we can generate a mass-maintaining set. Note that mass- 
maintaining sets are also self-maintaining, Gmm(Gsm(S)) = Gmm(S), which 
is a useful property, because Gsm(S) is easier to compute. As usual, the union 
□mm and intersection \1 MM of mass-maintaining sets S 1 ,S 2 are defined as 
S 1 U MM S 2 = Gmm^US^, S 1 \1 MM S 2 = G MM (Sir\S 2 ), respectively. Thus also 
the set of all mass-maintaining sets Omm forms a lattice (0mm, Umm, ^mm)- 
If S is mass-maintaining, its closure Gcl(S) is mass-maintaining, too (again, 
not valid for general reaction systems ). 

Finally, in consistent reaction systems , we can also generate uniquely an 
organization, (here, again, the largest organization that can be generated from 
a set) according to the following definition: 

Definition 11 Given a set of molecules C C J\A, we define G(C) as Gmm{Gcl{C)). 
We say that C generates the organization O = G(C). 

Equivalently G{C) = G mm{G sm{Gcl(C))) , which allows to compute the 
organization generated by a set more easily in three steps. Following the same 
scheme as before, the union U and intersection n of two organizations U and V 
is defined as the organization generated by their set-union and set-intersection: 
U U V = G(UUV), UnV = G(U (IV), respectively. Thus, for consistent 
reaction system , also the set of all organizations O forms a lattice (O, U, n). 
This important fact should be emphasized by the following lemma: 
Lemma 3 Given an algebraic chemistry (A4, 71) of a consistent reaction system 
and all its organizations O, then (O, U, n) is a lattice. 

Knowing that the semi-organizations and organizations form a lattice, and 
that we can uniquely generate an organization for every set, is a useful infor- 
mation. In order to find the whole set of organizations, it is impractical just 
to check all the possible sets of molecules. Instead, we can start by comput- 
ing the lattice of semi-organizations, and then test only those sets for mass- 
maintenance. Furthermore, if the semi-organizations form a lattice, we can 
start with small sets of molecules and generate their semi-organizations, while 
the Uso operator can lead us to the more complex semi-organizations. 
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As a summary, from a practical point of view, we calculate first the set 
of semi-organizations. If our system is a catalytic flow system , we automati- 
cally obtain the lattice of organizations. Otherwise we have to check for each 
semi-organization whether it is mass-maintaining or not. If we have a con- 
sistent reaction system , then we are assured to obtain a lattice, where there 
is a unique smallest and largest organization, and where we can easily obtain 
the intersection and union of two organizations from the graphical represen- 
tation of the lattice (see examples). For a general reaction system , the set of 
organizations does not necessarily form a lattice. Nevertheless, this set of orga- 
nizations represent the organizational structure of the reaction network, which 
can be visualized and which can provide a new view on the dynamics of the 
system by mapping the movement of the system in state space to a movement 
in the set of organizations, as will be shown in the following section. 

Dynamic Analysis 

The static theory deals with molecules M. and their reaction rules TZ, but not 
with the evolution of the system in time. To add dynamics to the theory, we 
have to formalize the dynamics of a system. In a very general approach, the 
dynamics is given by a state space X and a formal definition (mathematical 
or algorithmical) that describes all possible movements in X. Given an initial 
state x G X, the formal definition describes how the state changes over time. 
For simplicity, we assume a deterministic dynamical process, which can be 
formalized by a phase flow (X, (T t ) te ^) where (T t )t£R is a one-parametric group 
of transformations from X. T t (x ) denotes the state at time t of a system that 
has been in state x at t — 0. 

Connecting to the Static Theory A state x G X represents the state of 
a reaction vessel that contains molecules from M.. In the static part of the 
theory we consider just the set of molecular species present in the reaction 
vessel, but not their concentrations, spatial distributions, velocities, and so 
on. 

Now, given the state x of the reaction vessel, we need a function that maps 
uniquely this state to the set of molecules present. Vice versa, given a set of 
molecules A C M, we need to know, which states from X correspond to this 
set of molecules. For this reason we introduce a mapping called abstraction, 
from X to Ai , which maps a state of the system to the set of molecules that 
are present in the system being in that state. We require this mapping to be 
monotonically increasing on the number of molecules. In other words, if we 
have two states x x and x 2 , and the concentration of the molecule mi is higher 
in x 2 than in Xi; and if m x G <j>(x.i) then mi G 0(x 2 ). The exact mapping can 
be defined precisely later, depending on the state space, on the dynamics, and 
on the actual application. 
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The concept of instance is the opposite of the concept of abstraction. While 
</>(x) denotes the molecules represented by the state x, an instance x of a set 
A is a state where exactly the molecules from A are present according to the 
function <fi. 

Definition 12 We say that a state x G X is an instance of A C M., iff 0(x) = A 
In particular, we can define an instance of an organization O (if <f>(x) = O) 
and an instance of a generator of O (if G((f>(x)) = O). Loosely speaking we 
can say that x generates organization O. Note that a state x of a consistent 
reaction system , reactive flow system , and catalytic flow system is always an 
instance of a generator of one and only one organization O. This leads to the 
important observation that a lattice of organizations partitions the state space 
X, where a partition Xq implied by organization O is defined as the set of all 
instance of all generators of O: Xq = {x G X|G(0(x)) = O}. Note that as 
the system state evolves over time, the organization G(0(x(t))) generated by 
x(t) might change (see below, Fig. 2 and 3). 

Fixed Points are Instances of Organizations. Now we describe a the- 
orem that relates fixed points to organizations, and by doing so, underlining 
the relevancy of organizations. We will show that, given an ODE of a form 
that is commonly used to describe the dynamics of reaction systems, every 
fixed point of this ODE is an instance of an organization. We therefore as- 
sume in this section that x is a concentration vector x = (xi,x 2 , ■ ■ ■ ,x\m\), 
X = M} M \,Xi > where Xi denotes the concentration of molecular species i 
in the reaction vessel, and M. is finite. The dynamics is given by an ODE 
of the form x = Mv(x) where M is the stoichiometric matrix implied by the 
algebraic chemistry (M.,H) (reaction rules). v(x) = (i>i(x), . . . ,t> n (x)) G M' 7 ^ 
is a flux vector depending on the current concentration x, where \TZ\ denotes 
the number of reaction rules. A flux u, (x) describes the rate of a particular 
reaction j. For the function Vj we require only that Vj(x) is positive, if and 
only if the molecules on the left hand side of the reaction j are present in the 
state x, and otherwise it must be zero. Often it is also assumed that Vj(x) 
increases monotonously, but this is not required here. Given the dynamical 
system as x = Mv(x), we can define the abstraction of a state x formally 
by using a (small) threshold 6 > such that all fixed points have positive 
coordinates greater than 6. 

Definition 13 Given a dynamical system x = /(x) and let x be a state in X, 
then the abstraction 0(x) is defined by 

0(x) = {i\ Xl > Q,i G M}, <j):X^V{M), 9>0 (1) 

where Xi is the concentration of molecular species i in state x, and 6 is a 
threshold chosen such that it is smaller than any positive coordinate of any 
fixed point of x = /(x), xi > 0. 
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Setting G = is a safe choice, because in this case <p always meets the 
definition above. But for practical reasons, it makes often sense to apply a 
positive threshold greater zero, e.g., when we take into consideration that the 
number of molecules in a reaction vessel is finite. 

Theorem 1 Hypothesis: Let us consider a general reaction system whose re- 
action network is given by the algebraic chemistry (Ai,TZ) and whose dynamics 
is given by x = Mv(x) = /(x) as defined before. Let x' 6 X be a fixed point, 
that is, /(x') = 0, and let us consider a mapping (ft as given by Def. H21 which 
assigns a set of molecules to each state x. Thesis: 0(x') is an organization. 

Proof: We need to prove that 0(x') is closed and mass-maintaining: (a) 
Closure: Let us assume that 0(x') is not closed, then there exist a molecule k 
such that k ^ 0(x') and k is generated by molecules contained in 0(x'). Since, 
x' is a fixed point, /(x') = Mv(x') = 0. Now we decompose the stoichiometric 
matrix M into two matrices M + and M separating all positive from the 
negative coefficients, respectively, such that M = M + +M~ and M + v(x') > 0, 
and M^v(x') < (note that v(x') is always non-negative by definition). Let 
x k and x k be the k-th row of M + v'(x) and M~v'(x), respectively, which 
represent the inflow (production) and outflow (destruction) of molecules of 
type k. Note that x k + x k = x' k = (fixed point condition). Since we assumed 
that k is produced by molecules from 0(x'), x k must be positive, x k > and 
thus x k < 0. But this leads to a contradiction. In order to show this we 
have to differentiate the following two cases: (i) Assume that x' k = 0, then 
x k must be zero, too (by definition; and intuitively, because a molecule not 
present cannot vanish), (ii) Assume that < x' k < 9 (x' k > 6 needs not to be 
considered, because in that case k would be contained in </>(x).). Because we 
assumed that x' k is a coordinate of a fixed point, 6 must be smaller than x' k , 
which is again a contradiction. 

(b) Mass-maintaining: We have to show that 0(x') is mass-maintaining. 
Since x' is a fixed point Mv(x') = 0, which fulfills condition (3) of the defini- 
tion of mass-maintaining (Def. 6). From the requirements for the flux vector 
v, it follows directly that U(a-+b)(x') > for all A e Vm((P(x')), which ful- 
fills condition (1) of Def. 6. Following the same contradictory argument as 
before in (ii), x' k must be zero for k ^ 0(x), and therefore v^a^b)^') = for 
A ^ Vm^I'x)), which fulfills the remaining condition (2) of mass-maintaining, 
q.e.d. 

From this theorem it follows immediately that a fixed point is an instance 
of a closed set, a self-maintaining set, and of a semi-organization. Let us 
finally mention that even if each fixed point is an instance of an organization, 
an organization does not necessarily possess a fixed point. Further note that 
given an attractor A O X, there exists an organization O such that all points 
of A are instances of a generator of O. In fact, it might be natural to suppose 
that all points of an attractor are actually instances of O, yet it is not clear if 
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this is true for all systems or just for some. 
Movement from Organization to Organization 

ODEs and Movement in the Set of Organizations Not all system can 
be studied using ODEs. In particular a discrete system is usually only approx- 
imated by an ODE. In a discrete dynamical system, the molecular species that 
are present in the reaction vessel can change in time, e.g, as the last molecule 
of a certain type vanishes. In an ODE instead, this does not generally happen, 
where molecules can tend to zero as time tends to infinity. So, even if in reality 
a molecule disappears, in an ODE model it might still be present in a tiny 
quantity. The fact that every molecule ends up being present in (at least) 
a tiny quantity, generally precludes us to notice that the system is actually 
moving from a state where some molecules are present, to another state where 
a different set of molecules is present. Yet this is what happens in reality, 
and in this respect, an ODE is a poor approximation of reality. A common 
approach to overcome this problem is to introduce a concentration threshold 
0, below which a molecular species is considered not to be present. We use 
this threshold in order to define the abstraction 0, which just returns the set of 
molecules present in a certain state. Additionally, we might use the threshold 
to manipulate the numerical integration of an ODE by setting a concentra- 
tion to zero, when it falls below the threshold. In this constructive 
perturbation has to be greater than this threshold. 

Downward Movement Not all organizations are stable. The fact that 
there exits a flux vector, such that no molecule of that organization vanishes, 
does not imply that this flux vector can be realized when taking dynamics into 
account. As a result a molecular species can disappear. Each molecular species 
that disappears simplifies the system. Some molecules can be generated back. 
But eventually the system can move from a state that generates organization 
0\ into a state that generates organization 2 , with 2 always below 0\ 
(0 2 C Oi). We call this spontaneous movement a downward movement. 

Upward Movement Moving up to an organization above requires that a 
new molecular species appears in the system. This new molecular species can- 
not be produced by a reaction among present molecules (condition of closure). 
Thus moving to an organization above is more complicated then the movement 
down and requires a couple of specifications that describe how new molecular 
species enter the system. Here we assume that new molecular species appear 
by some sort of random perturbations or purposeful interference. We assume 
that a small quantity of molecules of that new molecular species (or a set 
of molecular species) suddenly appears. Often, in practice, the perturbation 
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(appearance of new molecular species) has a much slower time scale than the 
internal dynamics (e.g., chemical reaction kinetics) of the system. 

Example 2 Assume a system with two molecular species M. = {a, b} and the 
reactions 1Z = {a — > 2a, b — > 2b, a + b — > a, a — >,b — >}. All combinations of 
molecules are organizations, thus there are four organizations (Fig. 2). Assume 
further that the dynamics is governed by the ODE X\ = x% — x\, x 2 = x 2 — x\ — 
X\X 2 where X\ and x 2 denote the concentration of species a and b, respectively. 
We map a state to a set of molecules by using a small, positive threshold = 
0.1 > (Def. [T3J). Now, assume that the system is in state x = (0, 1) thus in 
organization {b}. If a small quantity of a appears (constructive perturbation), 
the amount of a will grow and b will tend to zero. The system will move 
upward to organization {a, b} in a transient phase, while finally moving down 
and converging to a fixed point in organization {a}. This movement can now 
be visualized in the lattice of organizations as shown in Fig. 2. 

Visualizing Possible Movements in the Set of Organizations In order 
to display potential movements in the lattice or set of organizations, we can 
draw links between organizations. As exemplified in Fig. 2 and Fig. 3, these 
links can indicate possible downward movements (down-link, blue) or upward 
movements (up- link, red). A neutral link (black line) denotes that neither the 
system can move spontaneously down, nor can a constructive perturbation 
move the system up. Whether the latter is true depends on the definition of 
"constructive perturbation" applied. For the example of Fig. 3 we defined a 
constructive perturbation as inserting a small quantity of one new molecular 
species. 

The dynamics in between organizations is more complex than this intuitive 
presentation might suggest, for example in some cases it is possible to move 
from one organization 0\ to an organization 2 , with 2 above (or below) 0\ 
without passing through the organizations in between 0\ and 2 . In Fig. 3 
this is the case for an upward movement from organization {a} to organization 
{a, b, c} caused by a constructive perturbation where a small quantity of c has 
been inserted. 

Organizations in Real Systems In preliminary studies JT] we have shown 
that artificial chemical reaction networks that are based on a structure-to- 
function mapping (e.g., ref. [TJ H2|) possess a more complex lattice of orga- 
nization than networks created randomly. From this observation we can al- 
ready expect that natural networks possess non-trivial organization structures. 
Our investigation of planetary photo-chemistries ^U] and bacterial metabolism 
(which will be published elsewhere), revealed lattices of organizations that van- 
ish when the networks are randomized, indicating a non-trivial structure. Here, 
in order to give an impression, we show a lattice of organizations obtained from 
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a model of the central sugar metabolism of E. coli by Puchalka and Kierzek jU] . 
The model consists of 92 species and 197 reactions, including gene expression, 
signal transduction, transport, and enzymatic activities. Figure 4 shows the 
lattice of organizations resulting from the original model, ignoring inhibiting 
interactions. The smallest organization, 0\, contains 76 molecules including 
the glucose metabolism and all input molecules. The input molecules, chosen 
according to ref. [D], include the external food set (Glcex, Glyex, Lacex) and 
all promoters (see suppl. material). Two other organizations, O3 and O4, con- 
tain the Lactose and Glycerol metabolism, respectively. Their union results in 
the largest organization O5 that contains all molecules. In summary, we can 
conclude, that the organizations found are biological meaningful, indicating a 
promising potential for future applications of this theory. 

Conclusion We presented a general way to define organizations in reaction 
systems, and proved a theorem that showed that fixed points are instances of 
those organizations. We also investigated the relative structure of the orga- 
nizations, opening up the door to visualize the dynamical movement of the 
system through organizations. Much work still needs to be done, and we be- 
lieve that this line of research is just in its infancy. Reaction systems with 
inhibiting interactions need to be investigated, as well as the structure inside 
organizations (i.e. the relation between attractive states and organizations). 
In fact the whole issue of movement among sets of molecules could only be 
introduced here. The concepts of attractive and stable organizations, which 
play a fundamental role in dynamical systems, has still to be formally defined. 
Yet we believe that our work represents a step forward towards a formal study 
of constructive dynamical system. 

We thank Florian Centler for providing the E. Coli. analysis. This work was supported by 
Federal Ministry of Education and Research (BMBF) Grant 0312704A to Friedrich-Schiller- 
University Jena. 
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FIGURE CAPTION 



Figure 1: Example with four species. Reaction network (a "2" means that 
two molecules are produced) and graphical representation of the lattice of all 
sets. The vertical position is determined by the set size. The solid lines depict 
the lattice of organizations. 



Figure 2: Example of an up- movement caused by a constructive perturba- 
tion, followed by a down- movement, (a) reaction network, (b) concentration 
vs. time plot of a trajectory, (c) lattice of organizations including trajectory. 



Figure 3: Lattice of organizations of Example 1, including up-links and 
down-links. Furthermore a trajectory is shown starting in organization {a, b, c, d} 
moving down to organization {a, b}. 



Figure 4: Lattice of organizations of a model of the central sugar metabolism 
of E.coli. jH]- In an organization, only names of new molecular species are 
printed that are not present in an organization below. 



16 




Figure 1: Example with four species. Reaction network (a "2" means that two 
molecules are produced) and graphical representation of the lattice of all sets. 
The vertical position is determined by the set size. The solid lines depict the 
lattice of organizations. 
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Figure 2: Example of an up- movement caused by a constructive perturbation, 
followed by a down- movement, (a) reaction network, (b) concentration vs. 
time plot of a trajectory, (c) lattice of organizations including trajectory. 
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Figure 3: Lattice of organizations of Example 1, including up-links and down- 
links. Furthermore a trajectory is shown starting in organization {a, b, c, d} 
moving down to organization {a, b). 
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organization 05 




Alio Lac LacZYmRNAl 
LacZYmRNA LacZ LacY 



organization 02 



CrpmRNA Crp 


organization Ol 






Glcex Glyex Lacex Glc6P Glc 
GlkPromCya .... 
[76 species total] 



Figure 4: Lattice of organizations of a model of the central sugar metabolism 
of E.coli. In an organization, only names of new molecular species are 
printed that are not present in an organization below. 
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